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Abstract 

Applications of Riemannian quantum geometry to cosmology have 
had notable successes. In particular, the fundamental discreteness un- 
derlying quantum geometry has led to a natural resolution of the big 
bang singularity. However, the precise mathematical structure under- 
lying loop quantum cosmology and the sense in which it implements the 
full quantization program in a symmetry reduced model has not been 
made explicit. The purpose of this paper is to address these issues, 
thereby providing a hrmer mathematical and conceptual foundation to 
the subject. 



1 Introduction 



In cosmology, one generally freezes all but a hnite number of degrees of 
freedom by imposing spatial homogeneity (and sometimes also isotropy). 
Because of the resulting mathematical simplifications, the framework pro- 
vides a simple arena to test ideas and constructions introduced in the full 
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theory both at the classical and the quantum levels. Moreover, in the clas- 
sical regime, the symmetry reduction captures the large scale dynamics of 
the universe as a whole quite well. Therefore, in the quantum theory, it 
provides a useful test-bed for analyzing the important issues related to the 
fate of classical singularities. 

Over the last three years, ramifications of Riemannian quantum geom- 
etry to cosmology have been investigated systematically. First, already at 
the kinematic level it was found that, thanks to the fundamental discrete- 
ness of quantum geometry, the inverse scale factor — and hence also the 
curvature — remains bounded on the kinematical Hilbert space [1]. Second, 
while classical dynamics is described by differential equations, the quantum 
Hamiltonian constraint can be interpreted as providing a difference equa- 
tion for the 'evolution' of the quantum state [2]. Furthermore, all quantum 
states remain regular at the classical big-bang; one can 'evolve' right through 
the point at which classical physics stops [3]. Third, the Hamiltonian con- 
straint together with the requirement — called pre-classicality — that the 
universe be classical at late times severely restricts the quantum state and, 
in the simplest models, selects the state uniquely [4]. There are also phe- 
nomenological models which allow us to study simple effects of quantum 
geometry leading to a behavior qualitatively different from the classical one 
[5]. Finally, the qualitative features are robust [6] and extend also to more 
complicated cosmological models [7] . These results are quite surprising from 
the perspective of the 'standard' quantum cosmology which was developed 
in the framework of geometro dynamics and, together, they show that, once 
the quantum nature of geometry is appropriately incorporated, the physical 
predictions change qualitatively in the Planck era. 

In spite of these striking advances, the subject has remained incomplete 
in several respects. First, in the existing treatments, certain subtleties which 
turn out to have important ramifications were overlooked and the underly- 
ing mathematical structure was somewhat oversimplified. This sometimes 
led to the impression that some of the physically desirable but surprising 
results arose simply because of ad-hoc assumptions. Second, the essential 
reasons why loop quantum cosmology is so different from the 'standard' 
quantum cosmology have not been spelled out. Third, while it is clear that 
the key constructions and techniques used in loop quantum cosmology are 
inspired by those developed in the full theory based on quantum geometry 
[8, 9, 10, 11, 12, 13, 14, 15], the parallels and the differences between the full 
theory and the symmetry reduced models have not been discussed in detail. 
In this paper, we will address these issues, providing a sounder foundation 
for the striking results obtained so far. The paper also has a secondary, 
pedagogical goal: it will also provide an introduction to quantum geome- 
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try and loop quantum gravity for readers who are not familiar with these 
areas. Our discussion should significantly clarify the precise mathematical 
structure underlying loop quantum cosmology and its relation to the full 
theory as well as to geometrodynamical quantum cosmology. However, we 
will not address the most important and the most difficult of open issues: a 
systematic derivation of loop quantum cosmology from full quantum gravity. 

Results of quantum cosmology often provide important qualitative lessons 
for full quantum gravity. However, while looking for these lessons, it is impor- 
tant to remember that the symmetry reduced theory used here differs from 
the full theory in conceptually important ways. The most obvious difference 
is the reduction from a field theory to a mechanical system, which eliminates 
the potential ultra-violet and infra-red problems of the full theory. In this 
respect the reduced theory is much simpler. However, there are also two 
other differences — generally overlooked — which make it conceptually and 
technically more complicated, at least when one tries to directly apply the 
techniques developed for the full theory. First, the reduced theory is usually 
treated by gauge fixing and therefore fails to be diffeomorphism invariant. As 
a result, key simplifications that occur in the treatment of full quantum dy- 
namics [14] do not carry over. Therefore, in a certain well-defined sense, the 
non-perturbative dynamics acquires new ambiguities in the reduced theory! 
The second complication arises from the fact that spatial homogeneity intro- 
duces distant correlations. Consequently, at the kinematical level, quantum 
states defined by holonomies along with distinct edges and triad operators 
smeared on distinct 2-surfaces are no longer independent. We will see that 
both these features give rise to certain complications which are not shared 
by the full theory. 

The remainder of this paper is divided in to four sections. In the second, 
we discuss the phase space of isotropic, homogeneous cosmologies; in the 
third, we construct the quantum kinematic framework; in the fourth we 
impose the Hamiltonian constraint and discuss properties of its solutions and 
in the fifth we summarize the results and discuss some of their ramifications. 



2 Phase space 

For simplicity, we will restrict ourselves to spatially flat, homogeneous, isotro- 
pic cosmologies, so that the spatial isometry group S will be the Euclidean 
group. Then the 3-dimensional group T of translations (ensuring homo- 
geneity) acts simply and transitively on the 3-manifold M. Therefore, M is 
topologically R 3 . Through the Cartan-Killing form on the Lie-algebra of the 
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rotation group, the Lie algebra of translations acquires an equivalence class 
of positive definite metrics related by an overall constant. Let us fix a metric 
in this class and an action of the Euclidean group on M. This will endow 
M with a fiducial flat metric °q a b. Finally, let us fix a constant orthonormal 
triad °ef and a co-triad °u l a on M, compatible with "q^. 

Let us now turn to the gravitational phase space in the connection vari- 
ables. In the full theory, the phase space consists of pairs (A l a , Ef) of fields 
on a 3-manifold M, where A l a is an SU(2) connection and Ef a triplet of 
vector fields with density weight 1 [16]. (The density weighted orthonormal 
triad is given by jEf, where 7 is the Barbero-Immirzi parameter.) Now, a 
pair (A'g 1 , E'f 1 ) on M will be said to be symmetric if for every s £ S there 
exists a local gauge transformation g : M — > SU(2), such that 

(s* A', s*E') = (g- l A'g + g^dg, g^E'g). (1) 

As is usual in cosmology, we will fix the local diffeomorphism and gauge free- 
dom. To do so, note first that for every symmetric pair (A', E') (satisfying 
the Gauss and diffeomorphism constraints) there exists an unique equivalent 
pair (A, E) such that 

A = c%\, E = p\f°q V (2) 

where c and p are constants, carrying the only non-trivial information con- 
tained in the pair (A',E'), and the density weight of E has been absorbed 
in the determinant of the fiducial metric. (Our conventions are such that 
[Tj, t j\ = eijkT k , i.e., 2iTi = o~i, where dj are the Pauli matrices.) 

In terms of p, the physical orthonormal triad ef and its inverse e l a (both 
of zero density weight) are given by: 

/°q _i 1 

- °ef = (sgnp) | 7 p| 2 ° e f, and e\ = (sgnp) 2 %\ (3) 

where q = det q a b = \ det jEf \ , sgn stands for the sign function and 7 is the 
Barbero-Immirzi parameter. As in the full theory, the Barbero-Immirzi pa- 
rameter 7 and the determinant factors are necessary to convert the (density 
weighted) momenta Ef in to geometrical (unweighted) triads ef and co- 
triads e\. The sign function arises because the connection dynamics phase 
space contains triads with both orientations and, because we have fixed a 
fiducial triad °ef, the orientation of the physical triad ef changes with the 
sign of p. (As in the full theory, we also allow degenerate co-triads which 
now correspond to p = 0, for which the triad vanishes.) 

Denote by As and Lg rav the subspace of the gravitational configuration 
space A and of the gravitational phase space r grav defined by (2). Tangent 
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vectors S to Tg rav are of the form: 

5 = (SA, SE), with 5A\ = (6c) %l 5E? = (Sp) °ef- (4) 

Thus, As is 1-dimensional and Tg rav is 2-dimensional: we made a restriction 
to symmetric fields and solved and gauge-fixed the gauge and the diffeo- 
morphism constraints, thereby reducing the local, infinite number of gravi- 
tational degrees of freedom to just one. 

Because M is non-compact and our fields are spatially homogeneous, 
various integrals featuring in the Hamiltonian framework of the full theory 
diverge. This is in particular the case for the symplectic structure of the full 
theory: 

fW*i, = j M d3x {SiAi(x)S 2 Et(x) - 5 2 A i a (x)5 1 Ef(x)) . (5) 

However, the presence of spatial homogeneity enables us to bypass this prob- 
lem in a natural fashion: Fix a 'cell' V adapted to the fiducial triad and 
restrict all integrations to this cell. (For simplicity, we will assume that 
this cell is cubical with respect to °q a b-) Then the gravitational symplectic 
structure f2 grav on r grav is given by: 

tt g rav(<5i, &) = J d 3 x (S 1 Ai(x)S 2 Ef(x) - 6 2 Ai(x)6 1 E?(x)) . (6) 

Using the form (4) of the tangent vectors, the pull-back of to Fg rav reduces 
just to: 

where V Q is the volume of V with respect to the auxiliary metric °q a b- (Had 
M been compact, we could set V = M and V Q would then be the total volume 
of M with respect to °q a b-) Thus, we have specified the gravitational part of 
the reduced phase space. We will not need to specify matter fields explicitly 
but only note that, upon similar restriction to symmetric fields and fixing of 
gauge and diffeomorphism freedom, we are led to a finite dimensional phase 
space also for matter fields. 

In the passage from the full to the reduced theory, we introduced a 
fiducial metric °q a b. There is a freedom in rescaling this metric by a constant: 
°qab k 2o q a b. Under this rescaling the canonical variables c,p transform via 
c i — ^ k~ l c and p k~ 2 p. (This is analogous to the fact that the scale factor 
a = y/\p\ in geometrodynamics rescales by a constant under the change 
of the fiducial flat metric.) Since rescalings of the fiducial metric do not 



238 



Mathematical structure of loop quantum cosmology 



change physics, by themselves c and p do not have direct physical meaning. 
Therefore, it is convenient to introduce new variables: 

c = VoC and p = V 3 p (8) 

which are independent of the choice of the fiducial metric °q a b- I n terms of 
these, the symplectic structure is given by 

n £av = dcAdp; (9) 

it is now independent of the volume V Q of the cell V and makes no refer- 
ence to the fiducial metric. In the rest of the paper, we will work with 
this phase space description. Note that now the configuration variable c 
is dimensionless while the momentum variable p has dimensions (length) 2 . 
(While comparing results in the full theory, it is important to bear in mind 
that these dimensions are different from those of the gravitational connec- 
tion and the triad there.) In terms of p, the physical triad and co-triad are 
given by: 



i , i . , , , i 



i 



e» = (sgnp)| 7P |-2 (V *°e°), and e\ = (sgnp)| 7 p|2 (V 5 % l a ) (10) 



Finally, let us turn to constraints. Since the Gauss and the diffeomor- 
phism constraints are already satisfied, there is a single non-trivial Scalar /Ha- 
miltonian constraint (corresponding to a constant lapse): 

6 

--2 c 2 sgnp v^pj + 8vrG C matt er = . (11) 



3 Quantization: Kinematics 

3.1 Elementary variables 

Let us begin by singling out 'elementary functions' on the classical phase 
space which are to have unambiguous quantum analogs. In the full theory, 
the configuration variables are constructed from holonomies h e (A) associated 
with edges e and momentum variables, from E(S,f), momenta E smeared 
with test fields / on 2-surfaces [13, 17, 18, 15]. But now, because of homo- 
geneity and isotropy, we do not need all edges e and surfaces S. Symmetric 
connections A in As can be recovered knowing holonomies h e along edges e 
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which lie along straight lines in M. Similarly, it is now appropriate to smear 
triads only across squares (with respect to "qab)- 1 

The SU(2) holonomy along an edge e is given by: 

h e (A) := V exp f A = cos — + 2 sin ^ (e ao Lo{) t 1 (12) 

J e 2 

1 

where fi £ (— oo, oo) (and fiVo is the oriented length of the edge with 
respect to °q a b)- Therefore, the algebra generated by sums of products of 
matrix elements of these holonomies is just the algebra of almost periodic 
functions of c, a typical element of which can be written as: 

?(c) = E^e«^ (13) 
j 

where j runs over a finite number of integers (labelling edges), [ij G R and 
£j £ C. In the terminology used in the full theory, one can regard a finite 
number of edges as providing us with a graph (since, because of homogeneity, 
the edges need not actually meet in vertices now) and the function g(A) as 
a cylindrical function with respect to that graph. The vector space of these 
almost periodic functions is, then, the analog of the space Cyl of cylindrical 
functions on A in the full theory [9, 10, 11, 13, 18]. We will call it the space 
of cylindrical functions of symmetric connections and denote it by Cyl s . 

In the full theory, the momentum functions E(S, f) are obtained by 
smearing the 'electric fields' Ef with an su(2)-valued function /' on a 2- 
surface S. In the homogeneous case, it is natural to use constant test func- 
tions p and let S be squares tangential to the fiducial triad °ef. Then, we 
have: 

E(S, f)= f Kbhdx a dx b = p V ^A S j (14) 
Js 

where = ry a b c i? c * and where A$j equals the area of S as measured by 
°q a b, times an obvious orientation factor (which depends on fi). Thus, apart 
from a kinematic factor determined by the background metric, the momenta 
are given just by p. In terms of classical geometry, p is related to the physical 
volume of the elementary cell V via 

V = \p\i. (15) 

1 Indeed, we could just consider edges lying in a single straight line and a single square. 
We chose not to break the symmetry artificially and consider instead all lines and all 
squares. 
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Finally, the only non-vanishing Poisson bracket between these elementary 
functions is: 

{g (A),p} = ^^^)e^. (16) 

3 

Since the right side is again in Cylg, the space of elementary variables is 
closed under the Poisson bracket. Note that, in contrast with the full theory, 
now the smeared momenta E(S, /) commute with one another since they are 
all proportional to p because of homogeneity and isotropy. This implies that 
now the triad representation does exist. In fact it will be convenient to use 
it later on in this paper. 

3.2 Representation of the algebra of elementary variables 

To construct quantum kinematics, we seek a representation of this algebra 
of elementary variables. In the full theory, one can use the Gel'fand theory 
to first find a representation of the C* algebra Cyl of configuration variables 
and then represent the momentum operators on the resulting Hilbert space 
[8, 9, 13, 18]. In the symmetry reduced model, we can follow the same 
procedure. We will briefly discuss the abstract construction and then present 
the explicit Hilbert space and operators in a way that does not require prior 
knowledge of the general framework. 

Let us begin with the C* algebra Cyl 5 of almost periodic functions on 
As which is topologically R. The Gel'fand theory now guarantees that there 
is a compact Hausdorff space -RBohr, the algebra of all continuous functions 
on which is isomorphic with Cylg. -RBohr is called the Bohr compactification 
of the real line As, and As is densely embedded in it. The Gel'fand theory 
also implies that the Hilbert space is necessarily L 2 (R-q^j. , d/i) with respect 
to a regular Borel measure fi. Thus, the classical configuration space As is 
now extended to the quantum configuration space -Reohr- The extension is 
entirely analogous to the extension from the space A of smooth connections 
to the space A of generalized connections in the full theory [8, 9, 13, 15] and 
came about because, as in the full theory, our configuration variables are 
constructed from holonomies. In the terminology used in the full theory, el- 
ements c of RBohr are 'generalized symmetric connections'. In the full theory, 
A is equipped with a natural, faithful, 'induced' Haar measure, which enables 
one to construct the kinematic Hilbert space and a preferred representation 
of the algebra of holonomies and smeared momenta [9, 10, 11, 12, 13]. 2 Sim- 
ilarly, -Reohr is equipped with a natural faithful, 'Haar measure' which we 

2 Recently, this representation has been shown to be uniquely singled out by the re- 
quirement of diffeomorphism invariance [19, 20, 21]. 
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will denote by /i D . 3 

Let us now display all this structure more explicitly. The Hilbert space 
^grav = L> 2 (RBohr,diJ, ) can be made 'concrete' as follows. It is the Cauchy 
completion of the space Cyl^ of almost periodic functions of c with respect 
to the inner product: 

<eH-|e*-l-)=5 M (17) 

(Note that the right side is the Kronecker delta, not the Dirac distribution.) 
Thus, the almost periodic functions A/^(c) := e* MC / 2 constitute an orthonor- 
mal basis in Hg rav . Cylg is dense in Hg rav , and serves as a common domain 
for all elementary operators. The configuration variables act in the obvious 
fashion: For all g\ and gi in Cylg, we have: 

(9i92)(c) = gi(c)g 2 (c) (18) 

Finally, we represent the momentum operator via 

V = -i 1 ^ J c , whence, (pg)(c) = ^ Y^fapj] M H (19) 

j 

where g £ Cyl s is given by (13) and, following conventions of loop quantum 
cosmology, we have set £pj = SttGH. (Unfortunately, this convention is 
different from that used in much of quantum geometry where Gh is set 
equal to £pi.) 

As in the full theory, the configuration operators are bounded, whence 
their action can be extended to the full Hilbert space Hg rav , while the mo- 
mentum operators are unbounded but essentially self-adjoint. The basis 
vectors jV^ are normalized eigenstates of p. As in quantum mechanics, let 
us use the bra-ket notation and write A/^(c) = (c|/x). Then, 

P \n) = ^\n) =pM. (20) 

Using the relation V = \p\ 3 ^ 2 between p and physical volume of the cell V 
we have: 

3 

v» = (^) 2 4» =vm. (21) 

3 ^Bohr is a compact Abelian group and d/i is the Haar measure on it. In non-relativistic 
quantum mechanics, using -RBohr one can introduce a new representation of the standard 
Weyl algebra. It is inequivalent to the standard Schrodinger representation and naturally 
incorporates the idea that spatial geometry is discrete at a fundamental scale. Nonethe- 
less, it reproduces the predictions of standard Schrodinger quantum mechanics within its 
domain of validity. (For details, see [22]). There is a close parallel with the situation in 
quantum cosmology, where the role of the Schrodinger representation is played by 'stan- 
dard' quantum cosmology of geometrodynamics. 
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This provides us with a physical meaning of fi: apart from a fixed constant, 
|/u| 3//2 is the physical volume of the cell V in Planck units, when the universe is 
in the quantum state |/x). Thus, in particular, while the volume V Q of the cell 
V with respect to the fiducial metric °q a b may be 'large', its physical volume 
in the quantum state \fi = 1) is (^/6) 3 ^ 2 £p V This fact will be important in 
sections 3.3 and 4.1. 

Note that the construction of the Hilbert space and the representation 
of the algebra is entirely parallel to that in the full theory. In particular, 
Cyl s is analogous to Cyl; -Reohr is analogous to A; to the spin network 
states Af a j,i (labelled by a graph g whose edges are assigned half integers j 
and whose vertices are assigned intertwiners I[23, 24, 18]). g are the analogs 
of configuration operators defined by elements of Cyl and p is analogous to 
the triad operators. In the full theory, holonomy operators are well-defined 
but there is no operator representing the connection itself. Similarly, are 
well defined unitary operators on Hg rav but they fail to be continuous with 
respect to fi, whence there is no operator corresponding to c on Hg rav . Thus, 
to obtain operators corresponding to functions on the gravitational phase 
space Fg rav we have to first express them in terms of our elementary variables 
Mn and p and then promote those expressions to the quantum theory. Again, 
this is precisely the analog of the procedure followed in the full theory. 

There is, however, one important difference between the full and the re- 
duced theories: while eigenvalues of the momentum (and other geometric) 
operators in the full theory span only a discrete subset of the real line, now 
every real number is a permissible eigenvalue of p. This difference can be 
directly attributed to the high degree of symmetry. In the full theory, eigen- 
vectors are labelled by a pair (e, j) consisting of continuous label e (denoting 
an edge) and a discrete label j (denoting the 'spin' on that edge), and the 
eigenvalue is dictated by j. Because of homogeneity and isotropy, the pair 
(e, j) has now collapsed to a single continuous label /x. Note however that 
there is a weaker sense in which the spectrum is discrete: all eigenvectors 
are normalizable. Hence the Hilbert space can be expanded out as a direct 
sum — rather than a direct integral — of the 1-dimensional eigenspaces of p; 
i.e., the decomposition of identity on Hs is given by a (continuous) sum 

rather than an integral. Although weaker, this discreteness is nonetheless 
important both technically and conceptually. In the next sub-section, we 
present a key illustration. 

We will conclude with two remarks. 
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i) In the above discussion we worked with c,p rather than the original 
variables c,p to bring out the physical meaning of various objects more di- 
rectly. Had we used the tilde variables, our symplectic structure would have 
involved V Q and we would have had to fix V a prior to quantization. The 
Hilbert space and the representation of the configuration operators would 
have been the same for all choices of V . However, the representation of the 
momentum operators would have changed from one V Q sector to another: A 
change V Q i— > k 3 V Q would have implied p i-> k~ 2 p. The analogous transfor- 
mation is not unitarily implementable in Schrodinger quantum mechanics 
nor in full quantum gravity. However, somewhat surprisingly, it is unitarily 
implementable in the reduced model. 4 Therefore, quantum physics does not 
change with the change of V Q . Via untilded variables, we chose to work with 
an unitarily equivalent representation which does not refer to V at all. 

ii) For simplicity of presentation, in the above discussion we avoided 
details of the Bohr compactification and worked with its dense space Cyl s 
instead. In terms of the compactification, the situation can be summarized 
as follows. After Cauchy completion, each element of Hg rav is represented 
by a square-integrable function f(c) of generalized symmetric connections. 
By Gel'fand transform, every element g of Cyl s is represented by a function 
g(c) on Cylg and the configuration operators act via multiplication on the 
full Hilbert space: (5132) (c) = <7i(c)02(c). The momentum operator p is 
essentially self- adjoint on the domain consisting of the image of Cyl s under 
the Gel'fand transform. 



3.3 Triad operator 

In the reduced classical theory, curvature is simply a multiple of the inverse 
of the square of the scale factor a = y/\p\. Similarly, the matter Hamiltonian 
invariably involves a term corresponding to an inverse power of a. Therefore, 
we need to obtain an operator corresponding to the inverse scale factor, or 
the triad (with density weight zero) of (10). In the classical theory, the triad 
coefficient diverges at the big bang and a key question is whether quantum 
effects 'tame' the big bang sufficiently to make the triad operator (and hence 
the curvature and the matter Hamiltonian) well behaved there. 

4 The difference from Schrodinger quantum mechanics can be traced back to the fact 
that the eigenvectors \p) of the Schrodinger momentum operator satisfy {p\p'} = 8(p,p') 
while the eigenvectors of our p in the reduced model satisfy {fJ.\(J,') = , the Dirac delta 
distribution being replaced by the Kronecker delta. In full quantum gravity, one also has 
the Kronecker-delta normalization for the eigenvectors of triad (and other geometrical) 
operators. Now the difference arises because there the eigenvalues form a discrete subset 
of the real line while in the symmetry reduced model they span the full line. 
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Now, in non-relativistic quantum mechanics, the spectrum of the op- 
erator f is the positive half of the real line, equipped with the standard 
Lesbegue measure, whence the operator 1/f is a densely-defined, self-adjoint 
operator. By contrast, since p admits a normalized eigenvector = 0) with 
zero eigenvalue, the naive expression of the triad operator fails to be densely 
defined on Wg rav . One could circumvent this problem in the reduced model 
in an ad-hoc manner by just making up a definition for the action of the 
triad operator on = 0). But then the result would have to be consid- 
ered as an artifact of a procedure expressly invented for the model and one 
would not have any confidence in its implications for the big bang. Now, as 
one might expect, a similar problem arises also in the full theory. There, a 
mathematically successful strategy to define the required operators already 
exits [14]: One first re-expresses the desired, potentially 'problematic' phase 
space function gular function of elementary variables and the volume 

function and then replaces these by their well-defined quantum analogs. It is 
appropriate to use the same procedure also in quantum cosmology; not only 
is this a natural approach but it would also test the general strategy. As in 
the general theory, therefore, we will proceed in two steps. In the first, we 
note that, on the reduced phase space Pg rav , the triad coefficient sgnp \p\~z 
can be expressed as the Poisson bracket {c, V 1 / 3 } which can be replaced 
by ih times the commutator in quantum theory. However, a second step is 
necessary because there is no operator c on Wg rav corresponding to c: one 
has to re-express the Poisson bracket in terms of holonomies which do have 
unambiguous quantum analogs. 



Indeed, on Pg rav , we have: 



sgn(p) 



\p\ 8vrG 7 



-tr 



{^T l h,{hT\V^ , (23) 



where 



v,? 



hi := V exp | / °e?AiTjdt J = exp(cT;) = cos | + 2 sin | Ti (2 1) 



is the holonomy (of the connection A l a ) evaluated along an edge along the 

1/3 

elementary cell V (i.e., an edge parallel to the triad vector °ef of length V 
with respect to the fiducial metric °q a b), and where we have summed over i 
to avoid singling out a specific triad vector. 5 We can now pass to quantum 



5 Note that, because of the factors hi and h i 1 in this expression, the length of the edge 
is actually irrelevant. For further discussion, see the remark at the end of this sub-section. 
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theory by replacing the Poisson brackets by commutators. This yields the 
triad (coefficient) operator: 



sgn(p) 



4i 
12* 



74 



tr (j2r%[h-\V l s]j 

(C ~ 1 C C " 1 C N 

sin -V* cos cos -Vs sin - 
2 2 2 2, 



(25) 



Although this operator involves both configuration and momentum op- 
erators, it commutes with p, whence its eigenvectors are again \fi). The 
eigenvalues are given by: 



sgn(p) 



lM> = 4- {V& - Vl'_\) | M > . (26) 



^2 VV+1 



where is the eigenvalue of the volume operator (see (21)). Next, we note 
a key property of the triad operator: It is bounded above! The upper bound 
is obtained at the value fx = 1: 

\p\dL = J^&. (27) 

This is a striking result because p admits a normalized eigenvector with zero 
eigenvalue. Since in the classical theory the curvature is proportional to 
in quantum theory, it is bounded above by (l2/j)£p^ . Note that h is essential 
for the existence of this upper bound; as H tends to zero, the bound goes to 
infinity just as one would expect from classical considerations. This is rather 
reminiscent of the situation with the ground state energy of the hydrogen 
atom in non-relativistic quantum mechanics, E = — (m e e 4 /2)(l/h), which 
is bounded from below because h is non-zero. 



In light of this surprising result, let us re-examine the physical meaning 
of the quantization procedure. (Using homogeneity and isotropy, we can 
naturally convert volume scales in to length scales. From now on, we will 
do so freely.) In the classical Poisson bracket, we replaced the connection 
coefficient c by the holonomy along an edge of the elementary cell V because 
there is no operator on Wg rav corresponding to c. Since the cell has volume 

1/3 

V with respect to the fiducial metric °q a bi the edge has length V ■ While 
this length can be large, what is relevant is the physical length of this edge 
and we will now present two arguments showing that the physical length 
is of the order of a Planck length. The first uses states. Let us begin 
by noting that, being a function of the connection, (matrix elements of) 
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the holonomy itself determine a quantum state A/^=i(A) = e 1 ? . In this 
state, the physical volume of the cell V is not V Q but (7/6) 3 / 2 -£ P1 . Hence, 
the appropriate physical edge length is (7/6) 1//2 ^pi, and this is only of the 
order of the Planck length. 6 The same conclusion is reached by a second 
argument based on the operator hf. Since e l ^\fi) = + 1) for any fx, the 
holonomy operator changes the volume of the universe by 'attaching' edges of 
physical length (7/6) 1 / 2 (|/i+l| 1 / 2 — l/uj 1 / 2 ) £p\. 7 These arguments enable us to 
interpret the quantization procedure as follows. There is no direct operator 
analog of c; only holonomy operators are well-defined. The 'fundamental' 
triad operator (25) involves holonomies along Planck scale edges. In the 
classical limit, we can let the edge length go to zero and then this operator 
reduces to the classical triad, the Poisson bracket {c, V 1 ^ 3 }. 

Since the classical triad diverges at the big bang, it is perhaps not sur- 
prising that the 'regularization' introduced by quantum effects ushers-in the 
Planck scale. However, the mechanism by which this came about is new and 
conceptually important. For, we did not introduce a cut-off or a regulator; 
the classical expression (23) of the triad coefficient we began with is exact. 
Since we did not 'regulate' the classical expression, the issue of removing 
the regulator does not arise. Nonetheless, it is true that the quantization 
procedure is 'indirect'. However, this was necessary because the spectrum 
of the momentum operator p (or of the 'scale factor operator' corresponding 
to a) is discrete in the sense detailed in section 3.2. Had the Hilbert space 
T^grav been a direct integral of the eigenspaces of p — rather than a direct 
sum — the triad operator could then have been defined directly using the 
spectral decomposition of p and would have been unbounded above. 

Indeed, this is precisely what happens in geometro dynamics. There, p 
and c themselves are elementary variables and the Hilbert space is taken to be 
the space of square integrable functions of p (or, rather, of a ~ If) 1 / 2 )- Then, 
p has a genuinely continuous spectrum and its inverse is a self-adjoint opera- 
tor, defined in terms the spectral decomposition of p and is unbounded above. 
By contrast, in loop quantum gravity, quantization is based on holonomies 
—the Wilson lines of the gravitational connection. We carried this central 
idea to the symmetry reduced model. As a direct result, as in the full theory, 
we were led to a non-standard Hilbert space Hg rav . Furthermore, we were 
led to consider almost periodic functions of c — rather than c itself — as 'el- 
ementary' and an operator corresponding to c is not even defined on Hg rav . 
All eigenvectors of p are now normalizable, including the one with zero eigen- 

6 Black hole entropy calculations imply that we should set 7 = -^2- to recover the 
standard quantum field theory in curved space-times from quantum geometry. [26] . 

7 The square-root of fx features rather than /1 itself because p corresponds to the square 
of the scale factor a and we chose to denote its eigenvalues by (■y/j,/Q)£p 1 . 
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value. Hence, to define the triad operator, one simply can not repeat the 
procedure of geometrodynamics. We are led to use the alternate procedure 
followed above. Of course one could simply invent a regularization scheme 
just for this symmetry reduced model. A key feature of our procedure is 
that it was not so invented; it is the direct analog of the procedure followed 
to address the same issue in the full theory [14]. 

Finally, let us return to the expression of the quantum operator (25). 
Since the fact that it is bounded is surprising, it is important to verify that 
the final result has physically reasonable properties. The first obvious re- 
quirement is that, since the triad coefficient sgnp/\p\2 is a function only of 
p, the triad operator should commute with p. A priori there is no guarantee 
that this would be the case. Indeed, the expression (25) of the triad oper- 
ator involves c as well. However, as we saw, this condition is in fact met. 
A second non-trivial requirement comes from the fact that the triad coef- 
ficient and the momentum are algebraically related in the classical theory: 
p- (sgnp/l^l 1 / 2 ) 2 = 1. A key criterion of viability of the triad operator is that 
this relation should be respected in an appropriate sense. More precisely, we 
can tolerate violations of this condition on states only in the Planck regime; 
the equality must be satisfied to an excellent approximation on states with 
large \i (i.e., with large volume). Is this the case? We have: 



Thus, up to order 0(^ 2 ), the eigenvalue of the triad operator is precisely 
sgnp^/ where p^ is the eigenvalue of p (see (21)). On states repre- 
senting a large universe 3> 1), the classical algebraic relation between 
the triad coefficient and p is indeed preserved to an excellent approximation. 
Violations are significant only on the eigen-subspace of the volume operator 
with eigenvalues of the order of or less, i.e., in the fully quantum regime. 

Remark: In (23), we used holonomies along edges of the elementary cell 
V. While this choice is natural because the cell is needed for classical con- 
siderations in any case, one might imagine using, instead, edges of length 
|// |^o • Had this been done, the replacement of (23) would have again 
provided an exact expression of the triad coefficient in the classical the- 
ory. However, to meet the second criterion above, one would be forced to 
choose \fi \ ~ 1 and we would be back with the 'natural' choice made above. 
Nonetheless, since there is a quantization ambiguity, the numerical coeffi- 
cients in the final results (e.g., the precise value of the upper bound of the 
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triad operator spectrum) should not be attached direct physical significance. 
In particular, for lessons for the full theory, one should use only the qualita- 
tive features and orders of magnitudes. The numerical values can be arrived 
at only through a systematic reduction of the full quantum theory, where 
the precise value of |// |(~ 1) should emerge, e.g., as the lowest eigenvalue of 
an appropriate geometric operator. We will return to this issue in the next 
section. 



4 Quantum dynamics: The Hamiltonian constraint 

Since the curvature is bounded above on the entire kinematical Hilbert space 
^grav one might expect that the classical singularity at the big bang would 
be naturally resolved in the quantum theory. In this section we will show 
that this is indeed the case. 

4.1 The quantum Hamiltonian constraint 

In section 1, we reduced the Hamiltonian constraint to (11). However, we 
can not use this form of the constraint directly because it is cast in terms 
of the connection c itself rather than holonomies. One can 'regulate' it in 
terms of holonomies and then pass to quantum theory. However, to bring 
out the close similarity of the regularization procedure with the one followed 
in the full theory [14], we will obtain the same expression starting from the 
expression of the classical constraint in the full theory: 

C gr av := / d^xNe- 1 (e ijk Fi b E^E bk - 2(1 + 7 2 )i^^4) 
J v 

= [ d 3 xNe ljk F l ab e- 1 E a ^E hk (29) 

Jv 

where e := y^det-E) sgn(det-E'). We restricted the integral to our cell V 
(of volume V a with respect to °q a b) and, in the second step, exploited the 
fact that for spatially flat, homogeneous models the two terms in the full 
constraint are proportional to each other (one can also treat both terms as 
in the full theory without significant changes [25]). Because of homogeneity, 
we can assume that the lapse N is constant and, for definiteness, from now 
onwards we will set it to one. 

As a first step in constructing a Hamiltonian constraint operator we have 
to express the curvature components F l ah in terms of holonomies. We will use 
the procedure followed in the full theory [14] (or in lattice gauge theories). 



A. Ashtekar, M. Bojowald and J. Lewandowski 



249 



Consider a square a>ij in the i-j plane spanned by two of the triad vectors 

1/3 

°ef, each of whose sides has length fi V with respect to the fiducial metric 
°q a b- 8 Then, 'the ab component' of the curvature is given by 

Fi b n = %i%i ( h ^l + {c^ )\ (30) 

The holonomy in turn can be expressed as 

h^ = h^hf°\h^y\hf°^ (31) 
where, as before, holonomies along individual edges are given by 

h \^) :=cos ^ + 2sin ^%. (32) 

Next, let us consider the triad term e~ l E 0-3 ' E hk in the expression of 
the Hamiltonian constraint. Since the triad is allowed to become degenerate, 
there is a potential problem with the factor e -1 . In the reduced model, e 
vanishes only when the triad itself vanishes and hence the required term 
tijk e~ l E a i E bk can be expressed as a non-singular function of p and the 
fiducial triads. In the full theory, the situation is more complicated and 
such a direct approach is not available. There is nonetheless a procedure to 
handle this apparently singular function [14]: one expresses it as a Poisson 
bracket between holonomies and the volume function as in section 3.3 and 
then promotes the resulting expression to an operator. To gain insight in 
to this strategy, here we will follow the same procedure. Thus, let us begin 
with the identity on the symmetry reduced phase space Fg rav : 

e tjk T l e~ x E a iE hk = -2(8vr 7 G ^V^)' 1 e abc % k c h^ o) {h^ o) -\ V} (33) 

where is the holonomy along the edge parallel to the kth basis vector 

1 /3 

of length fi V with respect to °q a b- Note that, unlike the expression (30) 
for F^ b , (33) is exact, i.e. does not depend on the choice of 

Collecting terms, we can now express the gravitational part of the con- 
straint as: 

C grav = -4(8vr 7 3 ^G)- 1 £ e ^ fe tr(^ ) / i ^ ) ^ ) - 1 ^ ) - 1 4^ ) {4^ ) - 1 , V}) 

ijk 

+Q(c 3 /i ) (34) 

8 In a model with non-zero intrinsic curvature, those edges would not form a closed 
loop. The issue of how to deal with intrinsic curvature is discussed in [27]. 
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where, the term proportional to identity in the leading contribution to 
F % ah in (30) drops out because of the trace operation and where we used 
e abc o^jio^o^k _ y/^ e u'fc. Note that, in contrast to the situation with triads 
in section 3.3, now the dependence on [i Q does not drop out. However, one 
can take the limit fi Q — > 0. Using the explicit form of the holonomies h^°\ 
one can verify that the leading term in (34) has a well-defined limit which 
equals precisely the classical constraint. Thus, now /j, — or the length of the 
edge used while expressing F a b in terms of the holonomy around the square 
a ij — plays the role of a regulator. Because of the presence of the curvature 
term, there is no natural way to express the constraint exactly in terms of 
our elementary variables; a limiting procedure is essential. This faithfully 
mirrors the situation in the full theory: there, again, the curvature term is 
recovered by introducing small loops at vertices of graphs and the classical 
expression of the constraint is recovered only in the limit in which the loop 
shrinks to zero. 

Let us focus on the leading term in (34). As in the full theory, this term 
is manifestly finite and can be promoted to a quantum operator directly. 
The resulting regulated constraint is: 

eg = ^(rV^r 1 ^ e ^ fc tr(^ ) ^ ) ^ ) ^ 1 / l ^ ) - 1 ^ ) [^ ) - 1 ,y]) 

ijk 

= 96,( 7 3 ^4i)^ S m 2 ^cos 2 ^ (35) 

xfsin^cos^-cos^sin^a 
V 2 2 2 2 / 

Its action on the eigenstates of p is 

Cfev» =3( 7 Vo4ir 1 (W- W)(|^ + 4 / u )-2|/i) + |/i-4 / u )). (36) 
On physical states, this action must equal that of the matter Hamiltonian 

— 87rGC mat ter- 

Now, the limit [i — > of the classical expression (34) exists and equals 
the classical Hamiltonian constraint which, however, contains c 2 (see (11)). 
Consequently, the naive limit of the operator C^&l also contains c 2 . How- 
ever, since c 2 is not well-defined on Hg rav , now the limit as [i Q — > fails 
to exist. Thus, we can not remove the regulator in the quantum theory of 
the reduced model. In the full theory, by contrast, one can remove the reg- 
ulator and obtain a well-defined action on diffeomorphism invariant states 
[14]. This difference can be directly traced back to the assumption of homo- 
geneity. 9 In the full theory, there is nonetheless a quantization ambiguity 

9 In the full theory, one triangulates the manifold with tetrahedra of coordinate volume 
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associated with the choice of the j label used on the new edges introduced to 
define the operator corresponding to F a b [29]. That is, in the full theory, the 
quantization procedure involves a pair of labels (e, j) where e is a continuous 
label denoting the new edge and j is a discrete label denoting the spin on 
that edge. Diffeomorphism invariance ensures that the quantum constraint 
is insensitive to the choice of e but the dependence on j remains as a quan- 
tization ambiguity. In the reduced model, diffeomorphism invariance is lost 
and the pair (e, j) of the full theory collapses into a single continuous label fi 
denoting the length of the edge introduced to define F a b. The dependence on 
\i Q persists — there is again a quantization ambiguity but it is now labelled 
by a continuous label fi . Thus, comparison of the situation with that in 
the full theory suggests that we should not regard Cg^av as an approximate 
quantum constraint; it is more appropriate to think of the /x -dependence in 
(35) as a quantization ambiguity in the exact quantum constraint. This is 
the viewpoint adopted in loop quantum cosmology. 

If one works in the strict confines of the reduced model, there does not 
appear to exist a natural way of removing this ambiguity. In the full theory, 
on the other hand, one can fix the ambiguity by assigning the lowest non- 
trivial j value, j = 1/2, to each extra loop introduced to determine the 
operator analog of F a b. This procedure can be motivated by the following 
heuristics. In the classical theory, we could use a loop enclosing an arbitrarily 
small area in the a-b plane to determine F a b locally. In quantum geometry, on 
the other hand, the area operator (of an open surface) has a lowest eigenvalue 
o-o = (v / 3t)/4^pi [17, 30] suggesting that it is physically inappropriate to 
try to localize F a b on arbitrarily small surfaces. The best one could do is to 
consider a loop spanning an area a Q , consider the holonomy around the loop 
to determine the integral of F a b on a surface of area a , and then extract an 
effective, local F a b by setting the integral equal to a F a b- It appears natural 
to use the same physical considerations to remove the quantization ambiguity 
also in the reduced model. Then, we are led to set the area of the smallest 
square spanned by to a , i.e. to set (7£t )^pi = a o, or {i Q = \/3/4. Thus, 
while in the reduced model itself, area eigenvalues can assume arbitrarily 
small values, if we 'import' from the full theory the value of the smallest 

/4 and writes the integral C(N) := J d 3 xN e^kF^e^ 1 E aj E bk as the limit of a Riemann 
sum, C(N) = ]im Mo ->o J2 V&NeijkF* b e~ 1 E aj E bk , where the sum is over tetrahedra. If 
we now replace F by a holonomy around a square a of length fi , F ~ fj,~ 2 (h a — 1), 
and the triad term by a Poisson bracket, e~ 1 EE ~ \i~ 1 h{hT x , V}, and pass to quantum 
operators, we obtain C(N) ~ lim trft a ft[/i _1 , V]. The /x factors cancel out but, in the 
sum, the number of terms goes to infinity as /i — ► 0. However, the action of the operator 
on a state based on any graph is non-trivial only at the vertices of the graph whence 
only a finite number of terms in the sum survive and these have a well defined limit on 
diffeomorphism invariant states. In the reduced model, because of homogeneity, all terms 
in the sum contribute equally and hence the sum diverges. 
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non-zero area eigenvalue, we are naturally led to set fi = \/3/4. We will do 
so. 

To summarize, in loop quantum cosmology, we adopt the viewpoint that 
(35), with \i = \/3/4, is the 'fundamental' Hamiltonian constraint operator 
which 'correctly' incorporates the underlying discreteness of quantum geom- 
etry and the classical expression (11) is an approximation which is valid only 
in regimes where this discreteness can be ignored and the continuum picture 
is valid. We will justify this assertion in section 4.3. 

4.2 Physical states 

Let us now solve the quantum constraint and obtain physical states. For sim- 
plicity, we assume that the matter is only minimally coupled to gravity (i.e., 
there are no curvature couplings). As in general non-trivially constrained 
systems, one expects that the physical states would fail to be normalizable 
in the kinematical Hilbert space 7i s = Wg rav <8> W ma tter ( see > e -S-> [31 j 18]). 
However, as in the full theory, they do have a natural 'home'. We again have 
a triplet 

Cyl 5 C H s C CyV s 

of spaces and physical states will belong to Cylg, the algebraic dual of Cylg. 
Since elements of Cylg need not be normalizable, we will denote them by 
(The usual, normalizable bras will be denoted by (*|.) 

It is convenient to exploit the existence of a triad representation. Then, 
every element of Cyl* s can be expanded as 

(*| = £>(^)<Ml (37) 

where 4> denotes the matter field and are the (normalized) eigenbras of 
p. Note that the sum is over a continuous variable \i whence (^| need not 
be normalizable. Now, the constraint equation 

(*l (Cfe] + SttGC-^L) 1 = (38) 
turns into the equation 

(W " V^Mcj), ii + 4^ D ) - 2(V^ - V^Wh fi) (39) 

+ (V3 Mo - 4 ^o) = -^87rG 7 3 ^4l ^matter 

for the coefficients ip(<f),fi), where C^°^ ter (/J,) only acts on the matter fields 
(and depends on (i via metric components in the matter Hamiltonian). Note 
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that, even though fi is a continuous variable, the quantum constraint is a 
difference equation rather than a differential equation. Strictly, (39) just 
constrains the coefficients ^{ip^n) to ensure that is a physical state. 
However, since each (fi\ is an eigenbra of the volume operator, it tells us how 
the matter wave function is correlated with volume, i.e., geometry. Now, if 
one wishes, one can regard p as providing a heuristic 'notion of time', and 
then think of (39) as an evolution equation for the quantum state of matter 
with respect to this time. (Note that p goes from — oo to oo, negative 
values corresponding to triads which are oppositely oriented to the fiducial 
one. The classical big-bang corresponds to p = 0.) While this heuristic 
interpretation often provides physical intuition for (39) and its consequences, 
it is not essential for what follows; one can forego this interpretation entirely 
and regard (39) only as a constraint equation. 

What is the fate of the classical singularity? At the big bang, the scale 
factor goes to zero. Hence it corresponds to the state \fi = 0) in Hg rav . So, 
the key question is whether the quantum 'evolution' breaks down at fx = 0. 
Now, the discrete 'evolution equation' (39) is essentially the same as that 
considered in the first papers on isotropic loop quantum cosmology [2, 3] 
and that discussion implies that the quantum physics does not stop at the 
big-bang. 

For completeness, we now recall the main argument. The basic idea is to 
explore the key consequences of the difference equation (39) which determine 
what happens at the initial singularity. Starting at fx = —ANfi for some 
large positive N, and fixing ip(cf), —ANix ) and ip(<j), (—4N+A)/j 10 ), one can use 
the equation to determine the coefficients ip(4>, (— 47V + 4n)fio) for all n > 1, 
provided the coefficient of the highest order term in (39) continues to remain 
non-zero. Now, it is easy to verify that the coefficient vanishes if and only if 
n = N. Thus, the coefficient ?/>(</>, /x = 0) remains undetermined. In its place, 
we just obtain a consistency condition constraining the coefficients ip(4>, A* = 
—4) and ^(4>,fi = —8). Now, since ip((j),fi = 0) remains undetermined, at 
first sight, it may appear that we can not 'evolve' past the singularity, i.e. 
that the quantum evolution also breaks down at the big-bang. However, 
the main point is that this is not the case. For, the coefficient tp(4>, n = 
0) just decouples from the rest. This comes about because, as a detailed 
examination shows, the minimally coupled matter Hamiltonians annihilate 
if)(4>, /i) for /i = [1, 25] and V^ a = V_ Mo . Thus, unlike in the classical theory, 
evolution does not stop at the singularity; the difference equation (39) lets us 
'evolve' right through it. In this analysis, we started at = —4N^i because 
we wanted to test what happens if one encounters the singularity 'head on'. 
If one begins at a generic /i, the 'discrete evolution' determined by (39) just 
'jumps' over the classical singularity without encountering any subtleties. 
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To summarize, two factors were key to the resolution of the big bang 
singularity: i) as a direct consequence of quantum geometry, the Hamiltonian 
constraint is now a difference equation rather than a differential equation 
as in geometrodynamics; and ii) the coefficients in the difference equation 
are such that one can evolve unambiguously 'through' the singularity even 
though the coefficient tp(<p, fj, = 0) is undetermined. Both these features are 
robust: they are insensitive to factor ordering ambiguities and persist in 
more complicated cosmological models [6, 7]. 

Next, let us consider the space of solutions. An examination of the 
classical degrees of freedom suggests that the freedom in physical quantum 
states should correspond to two functions just of matter fields (f>. The space 
of solutions to the Hamiltonian constraint, on the other hand is much larger: 
there are as many solutions as there are functions ip (</>,/•*) on an interval 
[// — 4/x ,// + 4/j, ), where jj! is any fixed number. This suggests that a 
large number of these solutions may be redundant. Indeed, to complete the 
quantization procedure, one needs to introduce an appropriate inner product 
on the space of solutions to the Hamiltonian constraint. The physical Hilbert 
space is then spanned by just those solutions to the quantum constraint 
which have finite norm. In simple examples one generally finds that, while 
the space of solutions to all constraints can be very large, the requirement 
of finiteness of norm suffices to produce a Hilbert space of the physically 
expected size. 

For the reduced system considered here, we have a quantum mechanical 
system with a single constraint in quantum cosmology. Hence it should be 
possible to extract physical states using the group averaging technique of 
the 'refined algebraic quantization framework' [31, 18, 32]. However, this 
analysis is yet to be carried out explicitly and therefore we do not yet have 
a good control on how large the physical Hilbert space really is. This issue 
is being investigated. 

The Hamiltonian constraint equation differs markedly from the Wheeler- 
DeWitt equation of geometrodynamics in the Planck regime because it cru- 
cially exploits the discreteness underlying quantum geometry. But one might 
expect that in the continuum limit fi — > — which, from the quantum ge- 
ometry perspective, is physically fictitious but nonetheless mathematically 
interesting — the present quantum constraint equation would reduce to the 
Wheeler-DeWitt equation. We will conclude this sub-section by showing 
that this expectation is indeed correct in a precise sense. 

To facilitate this comparison, it is convenient to introduce some notation. 
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Let us set 



C$£W,P) ■= o4i MR^OM]" = 87rGC mattcr (p)^(0,p) , (40) 



Then, the Wheeler-DeWitt equation can be written as 
2 
3 

where the prime denotes derivative with respect to p. If we now further set 

tip) := \po\V*(j*p )h^ ~ V s( P -Po)htl)^P) > 
with p = 7£p 1 /x G /6, our quantum constraint (39) reduces to: 

C^ v t^{4>,P) ■= -^4iPo 2 (V ; (0,P + 4p o )-2^((/>,p) + ^(0,p-4p o )) 

= 8^GC' mattcr (p)^(0,p). (41) 

From now on we will consider only those 'wave functions' ip((f>,p) which are 
smooth (more precisely, C 4 ) in their p dependence. Then, it follows that 

C^OM = 4^(^P)+4iO(Po)f'(^P) (42) 

+ 4i o(^) <(<M + 4i o(^)y/(<M + 4i o(j|)W,P) 

Hence, in the limit p — ► (i.e., fi — ► 0), we have 

C^j&p) ^ C££W,p) (43) 

whence the discrete equation (41) reduces precisely to the Wheeler-DeWitt 
equation (40). Put differently, it has turned out that (41) is a well-defined 
discretization of (40). 

One can also ask a related but distinct question: Is there a sense in 
which solutions to the Wheeler-DeWitt equation are approximate solutions 
to the 'fundamental' discrete evolution equation? The answer is again in the 
affirmative. Let us restrict ourselves to the part of the p-line where p 3> p , 
i.e., where the quantum volume of the universe is very large compared to the 
Planck scale. Consider the restriction, to this region, of a smooth solution 
tf)((j),p) to (40) and assume that it is slowly varying at the Planck scale in 
the sense that tp/tp' ~ s, tp/ft' ~ s 2 , etc, with p D <C s < p. Then, ip(<f>,p) 
is an approximate solution to the 'fundamental' quantum constraint (41) in 
the sense that: 



i + o(4) + o(4) + o(4?) + o(^; 

g Z pZ pi pi 



C^lmp) (44) 
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Note that, in contrast to the discussion about the relation between the two 
equations, we can not take the limit p Q — ► because we are now interested in 
the discrete evolution. The solution to the Wheeler-DeWitt equation is an 
approximate solution to the fundamental equation only to the extent that 
terms of the order 0(p 2 /s 2 ), 0(p 2 /p 2 ), 0(p 2 s/p ? '), 0(p 2 s 2 /p 4 ) are negligi- 
ble. 

We will conclude with three remarks. 

1) We saw in section 4.1 that the \i — ► limit of the quantum constraint 
operator c|rai does not exist on H^ rav . Yet, in the above discussion of 
the 'continuum limit', we were able to take this limit. The resolution of 
this apparent paradox is that the limit is taken on a certain sub-space of 
Cyl*, consisting of smooth functions of p and none of these states belong 
to Wg rav . Indeed, since elements of Hg rav have to be normalizable with 
respect to the inner product (17), they can have support only on a countable 
number of points; they cannot even be continuous. In particular, solutions 
to the Wheeler-DeWitt equation can not lie in 7i s ; they belong only to the 
enlargement Cyl* of H s . 

2) There is a close mathematical similarity between quantum cosmology 
discussed here and the 'polymer particle' example discussed in [22]. In that 
example, following the loop quantum gravity program, a new representa- 
tion of the Weyl algebra is introduced for a point particle in non-relativistic 
quantum mechanics. In this representation, the Weyl operators are unitarily 
implemented but weak continuity, assumed in the Von Neumann uniqueness 
theorem, is violated for one of the two 1-parameter unitary groups. As a 
result (although the position operator exists) the momentum operator — the 
generator of infinitesimal space translations — fails to exist. This is meant 
to reflect the underlying discreteness of geometry. The 'fundamental' quan- 
tum evolution is given by a difference equation. But there is a precise sense 
in which the standard Schrbdinger evolution is recovered in the regime of 
validity of non-relativistic quantum mechanics. The 'fundamental' discrete 
evolution is analogous to the present 'fundamental' quantum constraint (41) 
while the Schrodinger equation is the analog of the Wheeler-DeWitt equa- 
tion (40). Therefore, details of the polymer particle analysis provide good 
intuition for the 'mechanism' that allows loop quantum cosmology to be very 
different from the standard one in the Planck regime and yet agree with it 
when the universe is large compared to the Planck scale. 

3) As mentioned in section 1, in this paper we do not address the difficult 
issue of systematically deriving quantum cosmology from full loop quantum 
gravity. Indeed, since Cyl 5 <f_ Cyl, at first it seems it would be difficult to 
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relate the two theories. However, note that the physical states of the sym- 
metry reduced model are elements of Cylg while those of the full theory are 
elements of Cyl*, and Cyl s is contained in Cyl*: elements of Cyl^ are those 
distributions on the full quantum configuration space A which are supported 
only on the subspace As of symmetric connections [33]. In particular, solu- 
tions to the quantum constraint discussed in this section do belong to Cyl*. 
Therefore, it should be possible to recover such states by first considering 
the full quantum theory and then carrying out a symmetry reduction. 

4.3 Classical limit 

In section 4.1, we found that the gravitational part of the Hamiltonian con- 
straint could not be introduced by a straightforward 'quantization' of the 
classical constraint (11) because there is no direct operator analog of c on 
^grav We then followed the strategy adopted in the full theory to arrive at 

the expression (35) of Cgrai- To ensure that this is a viable quantization, 
we need to show that (35) does reduce to (11) in the classical limit. In this 
sub-section, we will carry out this task. 

For this purpose we will use coherent states peaked at points (c , N/j, ) 
of the classical phase space where N » 1 (i.e. the volume of the universe 
is very large compared to the Planck volume) and c Q <C 1 (late times, when 

— 1/3 

the extrinsic curvature is small compared to the fiducial scale V ). At 
such large volumes and late times, one would expect quantum corrections 
to become negligible. The question then is whether the expectation value 
of the quantum constraint Cgrav in these coherent states equals the classical 
constraint (11) modulo negligible corrections. If so, Cgf&v would have the 
correct classical limit. 

To construct a coherent state, we also have to specify the width d of the 
Gaussian (i.e., 'tolerance' for quantum fluctuations of p). Now, since the 
quantum fluctuations in the volume of the universe must be much smaller 
than the volume itself, d <C Nfi a and since we also want the uncertainty in 
c to be small, we must have fi <C d. A coherent state of the desired type is 
then given by: 

rifio) (45) 

(More precisely, |^) is the 'shadow' on the regular lattice /j, = nfi Q of the 
coherent state in Cyl* uniquely selected by the triplet (c , Nfi Q ,d). For 



!*> = £ 



^2 K 



exp(-((„-JV)^) expH((n-tfK)f) 
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details, see [22], Section 4.) Our task is to compute the expectation value 



(46) 



of the constraint operator (41): 

= 3(7Vo4i) _1 ( W - W) + 4 /"o) - 21^) + |/i - 4 Mo )) • 

Let us first calculate the expectation value. Setting e := /x /eZ, we have: 
(*\C\V) = exp(-^e 2 ((n / -iV) 2 + (n-iV) 2 ))e i tK-«V°(n , /i |(7|n /[ i ) 

n,n' 

= 3( 7 3 /Xo4i) _1 ^ exp(-ie 2 ((n' - iV) 2 + (n - iV) 2 ))e i tK-«W 

n,n' 

x(V(n+i)n ~ V(n-i)n ) {n'n \(\(n + A)n ) - 2|n/x ) + |(n - 4)/x )) 
= 3( 7 Vo4i)~ 1 [e 2 *^ Yl ex P("^ 2 (( n + 4 - iV) 2 + (n - N) 2 )) 

n ~ 

X (VlK ~ V (n-l)n ) 

-2^exp(-e 2 (n - N) 2 ) (V (n+1)lio - V {n _^ o ) 

n 

+e -2ic ^ £ exp( _ 1 e 2^ n _ 4 _ ^2 + (n _ ^2^ 

n 

X (^(n+1)^ - V(n-l)iJ,o) 

To simplify this expression further, we note that all three sums in this ex- 
pression are of the same form and focus on the first: 

J>xp(-ie 2 ((n + 4-iV) 2 + (n - N) 2 )) (V (n+1)(lo - V (n _ 1)/lo ) 



= e _4e2 V e-^C"-^) 2 



(V(n-l)/i _ ^(n-3)/i ) 



where we have completed the square in the exponential and shifted the sum- 
mation index by 2. To compute this sum, as in [22], we use the Poisson 
resummation formula 



Y^ e -^(n-Nf f{n) = Y, j e-^y- N ^f{y)e 2 ™y n dy. 



(47) 
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This integral can be evaluated using the steepest descent approximation (see 
Appendix). One obtains: 



£ e-^~ N ? f(n) = ^ E /(^ + 7?) e-^ +M (1 + 0((iV e )- 2 )) 

n n 

= ^f(N)(l + 0(e^ 2 ^) + 0((Ner 2 )) (48) 



where, in the last step, we used the fact that, since e <C 1, terms with n/0 
are suppressed by the exponential. (Note that iVe » 1 because N[i » d, 
i.e., the permissible quantum fluctuation in the volume of the universe is 
much smaller than the volume of the universe at the phase space point 
under consideration.) 

Finally, we can collect terms to compute the expectation value (46). 
Using = (V5r/e)(l + 0(e~^^)), we have 



(C) s = 3( 7 3 ^4i)- 1 [e- 4e2 e 2ic ^(V-iK-V-3K) 

"2(% +1)Mo - V {N _^ o ) + e^ 2 e-^°(V [N+ ^ o - V {N+1)flo ) 
x(l+ 0( e - 2 A 2 ) + 0((iVe) -2 )) 

= ^(7Vo)~V^pi/6 [e-^e 2 ^°((N - 1)1 - (TV - 3)§) 

- 2((7V + 1)1 - (N - 1)§) + e - 4e2 e- 2iCo ^((iV + 3)1 - (AT + 1)1) 
x(l + 0(e^ 2 / £2 ) + 0((iVe)- 2 )) 

= | (7 Vo) _1 \fwA^Tb (e~ 4e ' 2 e 2ic °»° - 2 + e^V 2 ^"") 

x(l + 0(e- 7r2 / £2 ) + 0((iVe)- 2 ) + O^ 1 )) 
= -6 7 - 2 c 2 v / P(l + 0((iVe)- 2 ) + 0(iV~ 1 ) + 0(e 2 ) + 0(c 2 )) (49) 



where we have set P := g7/V^pi-^ an d used the fact that c Q <C 1 and /u Q ~ 1 

2/2 

(we also dropped corrections of order e _7r ' e since they are always dwarfed 
by those of order e 2 ). Thus, the expectation value equals the classical con- 
straint (11) up to small corrections of order c 2 , ^/(Pe) 2 , ipi/P and e 2 . 
(Note that each of them can dominate the other corrections depending on 
the values of the different parameters.) Hence, (35) is a viable quantization 
of the classical expression (11). 
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5 Discussion 



Let us begin with a brief summary of the main results. In section 2, we car- 
ried out a systematic symmetry reduction of the phase space of full general 
relativity in the connection variables. In the spatially flat model considered 
here, our connection coefficient c equals the only non-trivial (i.e. dynam- 
ical) component of the extrinsic curvature (modulo a factor of 7) and our 
conjugate momentum p equals the only non-trivial metric component (mod- 
ulo sgnp). Hence, our symmetry reduced Hamiltonian description is the 
same as that of geometrodynamics. By contrast, we saw in section 3 that 
quantum theories are dramatically different already at the kinematic level. 
In loop quantum cosmology, the Hilbert space Wg rav is spanned by almost 
period functions of c while in geometrodynamics it would be spanned by 
square-integrable functions of c. The intersection between the two Hilbert 
spaces is only the zero element! In loop quantum cosmology, the fundamen- 
tal operators are p and = exp(i / uc/2); unlike in geometrodynamics, there 
is no operator corresponding to c itself. Although p is unbounded and its 
spectrum consists of the entire real line, all its eigenvectors are normaliz- 
able and the Hilbert space is the direct sum of the 1-dimensional sub-spaces 
spanned by the eigenspaces. In geometrodynamics, on the other hand, none 
of the eigenvectors of p is normalizable; the Hilbert space is a direct inte- 
gral of its 'eigenspaces'. This marked difference is responsible for the fact 
that, while the triad operator (which encodes the inverse of the scale factor) 
is unbounded in geometrodynamics, it is bounded in loop quantum cosmol- 
ogy. Consequently, in the state corresponding to the classical singularity, 
the curvature is large, but it does not diverge in loop quantum cosmology. 

In section 4 we discussed the Hamiltonian constraint, i.e., quantum dy- 
namics. Because there is no direct operator analog of c, we had to introduce 
the constraint operator Cgrav by an indirect construction. Here, we fol- 
lowed the strategy used in the full theory [14], expressing curvature F\ of 
the gravitational connection A l a in terms of its holonomies around suitable 
loops. In the full theory, one can take the limit as the loop shrinks to zero 
and obtain a well-defined operator on diffeomorphism invariant states. The 
reduced model, by contrast, fails to be diffeomorphism invariant and the 
operator diverges on Wg rav in the limit. Therefore, to obtain a well-defined 
operator, we used loops enclosing an area a a , the smallest non-zero quan- 
tum of area in quantum geometry. The resulting operator can be regarded 
as a 'good' quantization of the classical constraint function because it has 
the correct classical limit. The resulting quantum constraint equation has 
novel and physically appealing properties. First, it is a difference -rather 
than differential — equation and thus provides a 'discrete time evolution'. 
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Second, the coefficients in this difference equation are such that the 'evo- 
lution' does not break down at the singularity; quantum physics does not 
stop at the big-bang! This occurs without fine tuning matter or making it 
violate energy conditions. Furthermore, while in consistent discrete models 
the singularity is often 'avoided' because discrete 'time steps' are such that 
one simply leaps over the point where the singularity is expected to occur 
[36], here, one can and does confront the singularity head on only to find 
that it has been resolved by the quantum 'evolution'. Furthermore, these 
features are robust [6, 7]. However, near the big-bang, the state is 'extremely 
quantum mechanical,' with large fluctuations. Thus, the classical space-time 
'dissolves' near the big-bang. In this regime, we can analyze the structure 
only in quantum mechanical terms; we can no longer use our classical in- 
tuition which is deeply rooted in space-times and small fluctuations around 
them. 

In the Planck regime, the predictions of loop quantum cosmology are 
thus markedly different from those of standard quantum cosmology based 
on geometrodynamics. The origin of this difference can be traced back to the 
fact that while loop quantum cosmology makes a crucial use of the funda- 
mental discreteness of quantum geometry, standard cosmology is based on a 
continuum picture. One would therefore expect that the difference between 
the two would become negligible in regimes in which the continuum picture 
is a good approximation. We established two results to show that this expec- 
tation is indeed correct. First, there is a precise sense in which the difference 
equation of loop quantum cosmology reduces to the Wheeler-DeWitt differ- 
ential equation in the continuum limit. Second, in the regime far removed 
from the Planck scale, solutions to the Wheeler-DeWitt equation solve the 
difference equation to an excellent accuracy. Thus, the quantum constraint 
of loop quantum cosmology modifies the Wheeler-DeWitt equation in a sub- 
tle manner: the modification is significant only in the Planck regime and yet 
manages to be 'just right' to provide a natural resolution of the big-bang 
singularity. 

Next, let us re-examine the early papers on loop quantum cosmology in 
terms of the precise mathematical framework developed in this paper. In 
the present terminology, in the previous discussion one effectively restricted 
oneself just to periodic functions exp(inc/2), rather than almost periodic 
functions exp(i//c/2) considered here. Thus the gravitational Hilbert space 

S P 

Hgmv considered there is the rather small, periodic sub-space of the present 
Hp aT . While this restriction did have heuristic motivation, it amounted to 
forcing c to be periodic. 10 From the geometrodynamical perspective, the 



In the symmetry reduction, one began with the general geometric theory of invariant 
connections and found that components of a homogeneous connection transform as scalars 
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extrinsic curvature was made periodic (with a very large period) and it was 
then not surprising that the eigenvalues of the scale factor (and hence also 
the volume) operator could only be discrete. However, a careful analysis 
shows that the restriction to periodic functions can not be justified: periodic 
functions fail to separate the symmetric connections. Thus, in the earlier 
treatments, the space of configuration variables was 'too small' already at 
the classical level and this led to an artificial reduction in the size of the 
quantum state space. 

In the analysis presented here, c is not periodic. As a consequence, the 
spectrum of the volume operator is the entire real line. Yet, there is dis- 
creteness in a more subtle sense: all eigenvectors of the volume operator are 
normalizable. This is a direct consequence of the fundamental premise of 
loop quantum gravity that the quantum Hilbert space carries well defined 
operators corresponding to holonomies and not connections themselves. In 
the full theory, this feature does make the spectra of geometric operators 
discrete and their eigenvectors normalizable. Because of the homogeneity 
assumption, however, the first of these features is lost in loop quantum cos- 
mology but the second does survive. The surprising and highly non-trivial 
fact is that this is sufficient for several of the main results of earlier papers 
to continue to hold: i) the inverse scale factor is still bounded from above; 
ii) the Hamiltonian constraint is again a difference equation; and, iii) the 
coefficients in this equation are such that the singularity is resolved in the 
quantum theory. Furthermore, the current analysis provided a systematic 
approach to verify that the constraint operator Cgrav has the correct classi- 
cal limit and made its relation to the Wheeler-DeWitt operator more precise 
and transparent. However, as in earlier papers, the issue of finding the in- 
ner product on physical states is yet to be analyzed in detail. While the 
group averaging procedure [31, 18, 32] provides a natural avenue, a detailed 
implementation of this program has only begun. The issue of whether the 
'pre-classicality' condition selects unique quantum states, thereby providing 
a natural solution to the issue of initial conditions can be addressed system- 
atically only after one has a better control on the physical Hilbert space. 
Finally, the discussion of section 4.3 not only shows that the classical Ein- 
stein's equation is recovered in loop quantum cosmology in an appropriate 
limit but it also provides a systematic approach to the problem of finding 
quantum corrections to Einstein's equations. 11 These corrections are now 



under gauge transformations [34]. The appropriate, polymer theory for real- valued scalar 
fields was developed only recently [35] and requires, as in the present paper, the Bohr 
compactification of the real line. 

There appears to be a rather general impression that Einstein's equations are not 
modified in loop quantum gravity. As our discussion of 4.3 shows, this is not the case. It is 
true that we simply promoted the classical Hamiltonian constraint function to an operator. 
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being worked out systematically. 

We conclude with a general observation. The way in which the big- 
bang singularity is resolved has potentially deep implications on questions 
about the origin of the universe. For instance, the question of whether the 
universe had a beginning at a finite time is now 'transcended'. At first, 
the answer seems to be 'no' in the sense that the quantum evolution does 
not stop at the big bang. However, since space-time geometry 'dissolves' 
near the big-bang, there is no longer a notion of time, or of 'before' or 
'after' in the familiar sense. Therefore, strictly, the question is no longer 
meaningful. The paradigm has changed and meaningful questions must now 
be phrased differently, without using notions tied to classical space-times. 
A similar shift of paradigm occurred already with the advent of general 
relativity. Before Einstein, philosophers argued that the universe could not 
have a finite beginning because, if it did, one could ask what there was 
before. However, this question pre-supposes that space-time is an eternal, 
passive arena and matter simply evolves in it. With general relativity, we 
learned that space and time are 'born with matter', whence the question of 
'what was there before' is no longer meaningful. Loop quantum cosmology 
brings about a further shift of paradigm, weeding out certain questions that 
seemed meaningful in classical general relativity and requiring that they be 
replaced by more refined questions, formulated in the context of quantum 
space-times. 
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However, because there is no direct operator analog of c, this 'quantization' is subtle and 
even on semi-classical (coherent) states, sharply peaked at classical configurations, the 
expectation value of the constraint operator equals the classical constraint function with 
small but very specific quantum corrections. We expect that the same procedure can be 
applied in the full theory to obtain quantum corrections to full Einstein's equations. 
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A Method of Steepest Descent 

In this appendix, we will show that the method of steepest descent can be 
used to evaluate the right side of the Poisson resummation formula (47): 



To begin with, we will assume that / is analytic and return to the cases of 
interest at the end. 

Note first that the Fourier integral can be written as 



where g(z) = f(Nz) is the combination of volume eigenvalues, and 

h n (z) := -e 2 N{z - l) 2 + 2mnz . 

For large N, the integral can be evaluated by the method of steepest descent. 
For this, we have to first find the saddle points. In our case, this amounts to 
finding solutions z$ of h'(zo) = —2e 2 N(zQ — l) + 2mn = 0. This equation has 
a single solution, zq = 1 + mn/Ne 2 . Steepest paths through zq are defined 
by constant imaginary part Imh n (z) = Imh n (zo). If we write z = ^ + irj with 
real £, rj, we obtain —2e 2 N(£ — l)rj + 2irn^ = 2nn, which has the solutions 
£ = 1 or rj = irn/e 2 N. We can now deform the original integration along 
the real line to an integration along the steepest path rj = irn/e 2 N, i.e. 
z = t + iirn/e 2 N with t real. Next, let us change the integration variable to 
u by t = 1 + u/e\fN such that 



So far, everything was exact. In order to be able to compute the integral, 
we now use the Taylor expansion: 



g(t(u) +iTrn/e 2 N) = g(l + u/eVN + mn/e 2 N) = g(l + mn/e 2 N) 

+ug'(l +i7rn/e 2 N)/eVN + u 2 g"(l + iitn / e 2 N) / e 2 N 
+0(n 3 /e 3 iV5) 




n 



n 




u 2 = h n {z ) - K(z) = e 2 N{t - I) 2 . 



With dz = du/e\f~N we obtain 
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Then, the integral can be evaluated approximately as [37]: 

I g(z) e Nhn{z Uz = e^ivH g(l + iitn/e 2 N) f°° e N{h n {z Q )-v?) du 

J L J —oo 

/°° 
-oo 

/oo 
-oo 

= v^e-iiV-yi + mn/e 2 N)e~^ n2 l e2+2 ^ inN 
x(l+0({Ne)- 2 )), 

where the corrections of order (iVe)~ 2 are small because d <C N/j, a . (We 
assumed g"(zo) to be of the same order as g{zo), which is the case for the 
functions we are interested in here.) 

We can now return to the Poisson re-summation formula and compute 
the sum: 

Y^ e -e*{n-N? /(n) = J2 f e~^y- N ff(y)e 2 ™y n dy 

n n •* 

= [ f(Nz)e Nhn ^dz 

n 

x(l + 0((Ne)- 2 )) 
= v^^V(^)(l + 0(e-" 2 / £2 ) + 0((Ney 2 )) 

where, in the last step, we used the fact that terms with are suppressed 

by the exponential because t «1. 

Finally, let us address the fact that the functions f(y) of interest in the 
main text come from the eigenvalues of the volume operator which are of the 
type V nilo = (~f\n\n /6) 3 ' 2 and, because of the absolute value involved, f(y) 
fail to be analytic in y at y = 0. One can circumvent this potential problem 
by replacing them with functions f(y) which are the analytic continuations 
(to the upper half complex y-plane) of the restriction of f(y) to the positive 
real y-axis and then carry out the above calculation for / which, by construc- 
tion, are analytic. (Thus, for example, |y| 3 ^ 2 is replaced by (y 2 ) 3 ^ 4 .) The 
error, J2 n ex P(~ e2 ( n ~~ -^O 2 ) [f( n ) ~ /( n )L can be shown to be of the order 
0{e-( Ne ?) which is negligible compared to the corrections derived above. 
The error is so small because the function integrated on the right side of 
the Poisson resummation formula is strongly peaked at a very large positive 
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value y = N of y and the contribution from the negative half of the real y 
axis is extremely small. 
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